Resistance of superconducting nanowires connected to normal metal leads 
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We study experimentally the low temperature resistance of superconducting nanowires connected 
to normal metal reservoirs. We find that a substantial fraction of the nanowires is resistive, down to 
the lowest temperature measured, indicative of an intrinsic boundary resistance due to the Andreev- 
conversion of normal current to supercurrent. The results are successfully analyzed in terms of the 
kinetic equations for diffusive superconductors. 
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Superconducting nanowires are resistive down to very 
low temperatures due to dynamical changes of the macro- 
scopic phase (phase-slip). Thermally activated phase-slip 
(TAP) becomes more likely for reduced cross-sectional 
dimensions because the free-energy barrier scales with 
the area of the wire. Upon approaching T ^ 0, phase 
slip due to thermal activation disappears and resistivity 
persists only by macroscopic quantum tunneling through 
the free-energy barrier. These processes have recently 
been studied in suspended carbon nanotubes coated 
with a thin layer of a superconducting molybdenum- 
germanium (MoGe-) alloy/ and receive increased the- 
oretical attention. In a separate experiment^ ropes of 
single-walled carbon nanotubes show signs of supercon- 
ductivity, which should be strongly influenced by phase- 
slips as well. 

A second potential cause of low temperature resistance 
in superconducting wires is the penetration of a static 
electric field at normal-metal-superconductor interfaces. 
It reflects the conversion of current carried by normal 
electrons into one carried by Cooper-pairs via Andreev- 
reflection. It has been studied extensively close to the 
critical temperature, where it is related to quasiparti- 
cle charge imbalance. Although hardly experimentally 
studied, at very low temperatures a similar resistive con- 
tribution is expected to be present, reflecting a length 
of the order of the coherence length. Since the coher- 
ence length is a sizable portion of the resistive length of 
the nanowires it may contribute significantly to the mea- 
sured two-point resistance. Here we report experimental 
results on the resistance of narrow superconducting wires 
connected to normal metal leads(for short NSN). We find 
a strong contribution to the resistance due to the conver- 
sion processes, which is analyzed and understood in terms 
of the non-equilibrium theory for dirty superconductors. 

We chose to study samples (Fig^ made of supercon- 
ducting (S) aluminium (Al) because of its long coherence 
length. Our main interest is in the two-point resistance 
of the S-wire. Hence, we have chosen to work with thick 
and wide normal (N) contacts with a negligible contribu- 
tion to the normal state resistance. To minimize interface 
resistances due to electronic mismatch of both materials, 
we have chosen to work for N with bilayers of aluminium 
covered with thick normal metal (Cu). In such a ge- 




FIG. 1: SEM picture of a device (slightly misaligned), showing 
the coverage of the thin aluminium film with the thick Cu 
layer (except for one the measured devices are carefully lined 
up). The inset shows a schematic picture of an ideal device. 



ometry the superconducting aluminium wire is directly 
connected to normal aluminium. 

The S-wire of 100 nm thick Al is made by evapo- 
rating 99.999% purity Al at a rate of ~ 1 nm/s in a 
vacuum of 1 x 10~^ mbar during evaporation. Films 
made in the same way, have a residual resistance ra- 
tio, RRR = R300K /R4.2K of 7.5 indicative of the level 
of impurity scattering. Taking the phonon resistivity of 
pph{Al) = 2.7 /il^cm, the impurity resistivity is po = 
0.4 pQcm. Using ctq = N{0)e'^D and renormalized free- 
electron parameters: N{0) = 2.2 x 10^^ J~^m~^ and 
vf = 1.3 X 10^ m/s (Ref. 8), we find the elastic mean 
free path i of 100 nm, presumably limited by the thick- 
ness. The superconducting transition temperature of the 
100 nm film is 1.26 K, the usually enhanced value for 
aluminium thin films. 

The sample is made in one evaporation run using 
shadow evaporation. A 300 nm PMMA/ 500 nm PMMA- 
MAA lift-off mask defines the width and the length of 
the S-wire. The width is approximately 250 nm and 
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the length is varied from 1 /im to 4 /im. The wire and 
the aluminium of the normal reservoirs is deposited in 
one run, both 100 nm thick. The reservoirs are made 
normal by a 2"^ evaporation of 470 nm copper(Cu) on 
top. The Cu is evaporated at a rate of ~ 2.5 nm/s at 
a pressure of 1 x 10~^ mbar. The material-parameters 
are po = 0.4 /il^cm, £ = 165 nm, D = 66 x 10~^ m^/s 
using ^"(0) = 1.5 X 10"^^ J-^ir-^, and vf = 1.2 x 10^ 
m/s. The non-superconducting properties of the final 
reservoirs are confirmed by measuring the resistance of a 
100 nm Al/470 nm Cu bilayer down to the lowest temper- 
ature measured: 600 mK. No sign of superconductivity 
is observed. This is in agreement with the analysis of 
Tc for a bilayer using the Usadel equationA Only for a 
transparency of the Al/Cu interface lower than 0.1 the 
Tc would reach values above 600 mK. 0.1 is an unrealisti- 
cally low transparency for an in vacuum prepared Al/Cu 
interface. 

The two-point resistance is probed with a current of 1 
/I A (linear regime) in a ^ He system down to 600 mK. The 
voltage over the wire is measured with a lock-in amplifier 
at 133 Hz. 

We find that the normal state resistance of nominally 
identical wires scatters substantially (10% or more), pre- 
sumably due to grain-sizes compared to wire- widths. To 
allow a quantitative evaluation we have selected a set of 
wires with identical values of RRR = 3.2 ± 0.1. Fig- 
ure 121 shows the R — T measurements for this particular 
set of wires. Samples with different RRR values show 
qualitatively identical behavior. All wires show a finite 
remaining resistance as most striking result. 

Obviously, despite the differences in length, the resis- 
tances at low temperatures have identical values and fol- 
low the same trace. It indicates that the origin of this 
remaining resistance is likely due to the region in the S- 
wire next to the interface with the normal reservoir. The 
resistance at 600 mK is equal to a normal segment of the 
superconductor with a length of about 200 nm. 

Figure [2 (Inset) also shows that the critical temper- 
ature of the wire decreases linearly with increasing the 
inverse square of the wire length. 

The critical temperature of the wire should fol- 
low a straightforward Ginzburg-Landau analysis. For 
^gl(^) < ttL the bridge should become superconducting. 
This leads to an ~ 1/L^ onset-temperature according to: 
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For L ^ oo the normal contact can no longer depress Tc 
and we find the intrinsic critical temperature. For the 
studied wires it is found to be Tco = 1.26 K, identical to 
the independently determined values of the 100 nm Al 
film. 

From the normal state resistance we infer an impu- 
rity resistance of po = 1.1 pQcm^ in accordance with 
the RRR = 3.2. It is significantly higher than the 
po of the 100 nm Al film, leading to a lower diffusiv- 
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FIG. 2: (color online) Measured R — T curves for four dif- 
ferent bridge lengths. The intrinsic Tco is indicated by the 
vertical dashed line. The inset shows the measured critical 
temperature of the wire vs l/L^, which is used to determine 
Tco by letting L ^ oo. 



ity and elastic mean free path: D = 160 cm^/s, and 
^ = 37 nm, respectively. The resulting coherence length 
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124 nm. 
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Finally, we estimate a resistance contribution of 11 mQ 
due to the spreading resistance in the normal reservoirs 
at low temperatures, considerably less than the value we 
measure in Fig. [3 

Theoretically, since the studied nanowires show diffu- 
sive transport, the Usadel equations should apply to the 
system. It is most convenient to calculate the normal 
current for a given applied voltage difference, assuming 
linear response. Schmid and Schon^^ have shown that 
within this limit the normal current can be described 
with a variation, 5f{E^x)^ in the electronic distribution 
function f{E, x): 
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with X the coordinate along the length of the wire. 
6f{E^x) is determined from a Boltzmann-equation, 
which includes the conversion of electrons into Cooper- 
pairs, but ignores inelastic electron-phonon scattering 
(only relevant close to Tc): 
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The position-dependent spectral conductivity M{E^ x) 
consist of two parts Ni{E^x) = Re{G) and N2{E^x) = 
Re{F)^ with G the normal and F the anomalous Greens 
function: 



M{E,x) = iN,{E,x)f + {N2{E,x)y 



(4) 



Ni{E^x) is comparable to the standard BCS density of 
states. The applied voltage is taken into account via 
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FIG. 3: The calculated density-of- states 

Ni at various distances from the reservoirs 
(x = 100,200,300,400,500,1000 nm) for a 2 /xm long 
wire (t = OA,D = IGOcm'^/s and Ao = 192 /i eV). Note the 
exponentially small but finite sub-gap density-of-states in 
the middle (at x — l/xm; see inset). 
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FIG. 4: The voltage in the superconducting wire as a function 
of position for two different temperatures (t = 0.4 and 0.9). 
At t — 0.9 the wire behaves as a normal metal and for t = 0.4 
the voltage is clearly present to a depth ^ (wire length 2 /xm 
with D = 160 cm^/s and Ao = 192 /neV). The inset shows 
the position dependent normal currents and supercurrents 
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Hence, the normal metal leads are taken as equilibrium 
reservoirs. The strength of the pairing interaction A{x) 
in Eq. [2lis determined by solving the Usadel equation in 
the Matsubara representation: 
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in conjunction with the self-consistency equation 
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The so-called proximity angle, ^, parameterizes the nor- 
mal Green's function G = cos^ and the anomalous 
Green's function F = sinO. The spectral functions are 
calculated by again solving the Usadel equation but now 
for CO —iE. The large normal metal reservoirs impose 
the boundary condition 6{x = 0, L) = 

Our main interest is the question how the current con- 
version process contributes to the resistance. First of 
all, the presence of decaying normal electron states sup- 
presses the gap in the density of states. 

In Fig. 13 the density-of-states Ni{E) is given for sev- 
eral positions along the wire of L = 2 /im, D = 160 
cmVs, Ao = 192 /ieV, and T/T^ = 0.4. Clearly, mov- 
ing away form the normal leads the density of states re- 
sembles more and more the well-known BCS density of 



states. Note however, that a finite sub-gap value remains 
in the middle {x = 1/im) even for very long wires. This is 
an intrinsic result for any NSN system and it is not due 
to current-flow, since this result is calculated in thermal 
equilibrium. (The back-action of the current-flow on the 
spectral properties can be neglected in the linear response 
regime). 

In Fig. El we show the results of a calculation of 
the voltage as a function of position along the wire for 
two different temperatures: t = 0.4, and t = 0.9 with 
D = 160 cm^/s, and Aq = 192 fieV. At the tempera- 
ture close to the transition temperature, the electric field 
penetrates the sample completely and the resistance is 
close to the normal state value. At low temperatures, 
the electric field still penetrates the superconductor over 
a finite length, leaving a middle piece with hardly any 
voltage drop.^^ The penetration length is of the order of 
the coherence length. The inset shows the position de- 
pendent normal currents and supercurrents illustrating 
the current conversion processes. 

In Fig. O a comparison is made between the calculated 
resistance as a function of temperature and the measure- 
ment for a L = 2 jam wire. The calculation is done with 
parameters D = 160 cm^/s, as determined from the im- 
purity resistivity, and A(0) = 1.764A:bTc = 192 /ieV with 
Tc = 1.26 K determined from the length dependence. 
These parameters have been determined independently. 
Without any fitting parameter, the agreement between 
the model (dots) and the experiment (data points: open 
symbols) is encouragingly good. Only at the lower tem- 
peratures the observed resistance is slightly less than the 
theoretically predicted values. 

Apparently, the model is overestimating the remaining 
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FIG. 5: (color online) The measured R-T curve for the 2 /xm 
bridge together with the model calculation using the bound- 
ary condition 0{x = 0,L) = (dots), and using the bound- 
ary condition dxO{x = 0,L) = 0{x = 0,L)/a (triangles) with 
a = 75 nm. The value for T=0 is 0.3255 for the hard boundary 
conditions (dots), and 0.2609 for the soft boundary conditions 

resistance below T = 0.8 K . Since there is little freedom 
left, we have very few options to remedy this discrep- 
ancy. The most likely option is that the rigid boundary 
conditions imposed at the interfaces should be relaxed. 
There is a finite possibility for superconducting correla- 
tions to extend into the normal metal reservoirs, which 
would mean that the boundary condition 6{x = 0, L) = 
is too rigid. Since the correlations will extend into a 3- 
dimensional volume we assume that using the boundary 
conditions dxO{x = 0, L) = 0{x = 0,L)/a i.e. a decay 
over a fixed characteristic length a, is a realistic assump- 
tion. It assumes a geometric dilution of the correlations. 
The result is shown in Fig. [Slby filled triangles. The best 
agreement between measurement and calculation model 
is obtained for a = 75 nm. This value appears reasonable 
for a decay length since it is comparable to the dimen- 
sions (100 nm by 250 nm) of the wire which emits into the 
reservoirs. The shortest bridge shown in Fig. [21 is found 
to have a further reduction in resistance, which we at- 



tribute to a small misalignment as shown in Fig. ^ Note 
however that the model predicts that NSN devices will 
continue to be resistive down to T = K. For the rigid 
boundary conditions we find R/R4.2K = 0.3255 and for 
the soft boundary conditions we find R/R4.2K = 0.2609 
An early indication of the low temperature boundary 
reported here is given by Harding et al}^. They stud- 
ied the resistance of thick sandwiches of Pb-Cu-Pb. By 
subtracting the contribution to the resistance of the Cu- 
layers they identified a remaining boundary resistance 
which depended on the mean free path in the Pb lay- 
ers. In later work by Hsiang and Clarkei such a resis- 
tance appeared to be unobservable, in contrast however 
to more recent work by Gu et al.^^ The geometry of our 
sample, with a negligible contribution to the resistance 
from the reservoirs, allows us to measure only the re- 
sistance due to the conversion processes in the super- 
conductor. In response to the work of Harding et al. 
Krahenbiihl and Watts- Tobin^^ derived an analytical ex- 
pression for the effective length of the boundary resis- 
tance: xq = ^ ^7T^olf{T) with £ the mean free path for 

elastic scattering, ^0 the BCS coherence length and / a 
function of temperature. For T = K, the function / is 
of order 1. In contrast to our model Ref. [l^ allows for 
one-dimensional diffusion in N, rather than treating N as 
an equilibrium reservoir. 

In conclusion, we have shown that the resistance of 
a superconducting nano-wire connected to normal leads 
has a finite DC resistance down to very low temperatures. 
The microscopic theory describes the data very well and 
provide a detailed image of the conversion from normal 
current into super current, over a few coherence lengths. 
The results emphasize the important role of the length 
of the wires in relation to the nature of the contacts.^ It 
also explains results obtained with diffusion-cooled hot- 
electron bolometers in which normal leads are used to 
provide rapid out-diffusion of hot electrons from a super- 
conducting wire^^. 
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